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Abstract 

In this paper, we study the design of pulse sequences for NMR spectroscopy as a problem of time 
optimal control of the unitary propagator. Radio frequency pulses are used in coherent spectroscopy 
to implement a unitary transfer of state. Pulse sequences that accomplish a desired transfer should 
be as short as possible in order to minimize the effects of relaxation and to optimize the sensitivity of 
the experiments. Here, we give an analytical characterization of such time optimal pulse sequences 
applicable to coherence transfer experiments in multiple-spin systems. We have adopted a general 
mathematical formulation, and present many of our results in this setting, mindful of the fact that 
new structures in optimal pulse design are constantly arising. Moreover, the general proofs are 
no more difficult than the specific problems of current interest. From a general control theory 
perspective, the problems we want to study have the following character. Suppose we are given 
a controllable right invariant system on a compact Lie group, what is the minimum time required 
to steer the system from some initial point to a specified final point? In NMR spectroscopy and 
quantum computing, this translates to, what is the minimum time required to produce a unitary 
propagator? We also give an analytical characterization of maximum achievable transfer in a given 
time for the two-spin systems. 



1 Introduction 

Many areas of spectroscopic fields, such as nuclear magnetic resonance (NMR), electron magnetic 
resonance and optical spectroscopy rely on a limited set of control variables in order to create de- 
sired unitary transformations ^, In NMR, unitary transformations are used to manipulate 
an ensemble of nuclear spins, e.g. to transfer coherence between coupled spins in multidimensional 
NMR-experiments || or to implement quantum- logic gates in NMR quantum computers ||. How- 
ever, the design of a sequence of radio-frequency pulses that generate a desired unitary operator is 
not trivial |J . Such a pulse sequence should be as short as possible in order to minimize the effects 
of relaxation or decoherence that are always present. So far, no general approach was known to de- 
termine the minimum time for the implementation of a desired unitary transformation Here we 
give an analytical characterization of such time optimal pulse sequences related to coherence transfer 
experiments in multiple spin systems. We determine, for example, the best possible in-phase and 
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anti-phase || [[(J coherence transfer achievable in a given time. We show that the optimal in-phase 
transfer sequences improve the transfer efficiency relative to the isotropic mixing sequences fll|| and 
demonstrate the optimality of some previously known sequences. 

During the last decade the questions of controllability of quantum systems have generated con- 
siderable interest Jll| [ijj . In particular, coherence or polarization transfer in pulsed coherent spec- 
troscopy has received lot of attention || Algorithms for determining bounds quantifying the 
maximum possible efficiency of transfer between non-Hermitian operators have been determined j(| . 
There is utmost need for design strategies for pulse sequences that can achieve these bounds. From 
a control theory perspective, this is a constructive controllability problem |L4]]. At the same time it 
is desirable that the pulse sequences be as short as possible so as to minimize the relaxation effects. 
This naturally leads us to the problem of time optimal control, i.e. given that there exist controls 
that steer the system from a given initial to final state, we would like to determine controls that 
achieve the task in minimum possible time |T^, [l^] . 

In non-relativistic quantum mechanics, the time evolution of a quantum system is defined through 
the time-dependent Schrodinger equation 

U{t) = -iH(t)U(t), U(0)=I, 

where H(t) and U(t) are the Hamiltonian and the unitary displacement operators, respectively. In 
this paper, we will only be concerned with finite-dimensional quantum systems. In this case, we can 
choose a basis and think of H(t) as a Hcrmitian matrix. We can split the Hamiltonian 

m 

i=i 

where is the part of Hamiltonian that is internal to the system and we call it the drift or free 
Hamiltonian and ^<=i v i(t)Hi is the part of Hamiltonian that can be externally changed. It is called 
the control or rf Hamiltonian. The equation for U(t) dictates the evolution of the density matrix 
according to 

p(t) = U(t)p(0)U^(t). 

The problem we are ultimately interested in is to find the minimum time required to transfer the 
density matrix from the initial state po to a final state pf- Thus, we will be interested in computing 
the minimum time required to steer the system 

m 

U = -HH i + Y i v i H i )U, (1) 

from identity, U(Q) = /, to a final propagator Uf- 

In the following section we establish a framework for studying such problems. For reasons 
suggested before our approach is more general than the current application requires, but this added 
generality does not complicate the development. 



2 Preliminaries 



We will assume that the reader is familiar with the basic facts about Lie groups and homogeneous 
spaces Q. Throughout this paper, G will denote a compact semi-simple Lie group and e its identity 
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element (we use I to denote the identity matrix when working with the matrix representation of the 
group). As is well known there is a naturally defined bi- invariant metric on G, given by the Killing 
form. We denote this bi-invariant metric by <, ><?. Let K be a compact closed subgroup of G. We 
will denote by L(G) the Lie algebra of right invariant vector fields on G and similarly L(K) the Lie 
algebra of right invariant vector fields on K . There is a one to one correspondence between these 
vector fields and the tangent spaces T e (G) and T e (K), which we denote by g and t respectively. 
Consider the direct sum decomposition g = m + 1 such that m = t 1 - with respect to the metric. 

To fix ideas, let G = SU (n) and g — su(n) be its associated Lie algebra ofnxn traceless skew- 
Hermitian matrices. Then < A, B >g= tr(A> B), A,B£ su(n) (which is proportional to the Killing 
metric) represents a bi-invariant metric on SU(n). 

It is well known that the (right) coset space G/K = {KU : U 6 G} (homogeneous space) admits 
the structure of a differentiable manifold Let 7r : G — > G/K denote the natural projection map. 
Define o E G/K by o — 7r(e). Given the decomposition g = m + 6, there exists a neighborhood 
of € m which is mapped homeomorphically onto a neighborhood of the origin o 6 G/K by the 
mapping 7r o exp| m . The tangent space plane T Q (G/K) can be then identified with the vector 
subspace m. The geometry of homogeneous space will play an essential part in determining the 
shortest possible times for transfers. 

The Lie group G acts on its Lie algebra g by conjugation Ado '■ g — > g (called the adjoint action). 
This is defined as follows. Given U G G, leg, then 

dt 

Once again to fix ideas if G = SU(n) and U £ G, A 6 su(n), then Adu(A) = U^AU. We use the 
notation 

Ad K (X) = |J Ad k (X). 

k£K 

If the homogeneous space G/K is a Riemannian symmetric space ||, the Lie algebra decom- 
position g = m + 6 (sec for properties of these orthogonal involutive Lie algebras) satisfies the 
commutation relation 

[6, t] C 6, [m, t] c m, [m, m] C t. 

If f) is a subalgebra of g contained in m, then f) is abelian because [m, m] € g. It is well known || 
that: 

Theorem 1 If f) and \)' are two maximal abelian subalgebras of m, then 

1. There exists an element £ € f) whose centralizer in m is just f). 

2. There is an element k 6 K such that Adk(i)) = fj'- 
3- m = [j keK Ad k (l)). 



Thus the maximal abelian subalgebras of m are all Adx conjugate and in particular they have 
the same dimension. The dimension will be called the rank of the symmetric space G/K and the 
maximal abelian subalgebras of m are called the Cartan subalgebras of the pair (g,t). We will see 
in what follows that the structure of the time optimal control depends on the rank in an important 
way. We state a useful corollary of the above the theorem ||] . 



Corollary 1 Let G/K be a Riemannian symmetric space. Let t) be a Cartan subalgebra of the pair 
(g, t) and define A = exp(f)) C G. Then G = KAK. 
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Proof: G = KP, where P = exp(m) = exp([J keK Ad k (i))) = \JkeK Ad k (exp{t))) = \J keK Ad k (A) C 
KAK. Now G = KKAK = KAK. Q.E.D. 

Note the space G/K is a union of maximal abelian subgroups Ad k (A), called maximal tori. 
Assumption 1 Let U € G and let the control system 

m 

17 = [X d + ViXi]U, 17(0) = e (2) 

be given. Please note we are working with the matrix representation of the group. We use 
{Xd, Xi, . . . , X, n } la to denote the Lie algebra generated by {Xd, X\, . . . , X m }. We will assume 
that {Xd, X\, . . . , X m }LA = 0, and since G is compact, it follows that the system (0) is controllable 
j|. Let t = {Xi}i,A and K — exp{Xi}LA be the closed compact group generated by {Xi}. Given 
the direct sum decomposition q = m + t where m = t - 1 with respect to the bi- invariant metric <, >g, 
let Xd G m. We will assume that Ada (ta) C m, in which case one says the homogeneous space G/K 
is reductive. All our examples will fall into this category. 

Notation: Let C denote the class of all locally bounded measurable functions defined on the interval 
[0, oo) and taking value in R m . C[0, T] denotes their restriction on the interval [0, T]. We will assume 
throughout that in equation (||), v — (v\,V2, ■ ■ ■ ,v m ) G C. Given v G C, we use U(t) to denote the 
solution of equation such that U(0) = e. If, for some time t > 0, U(t) = U' , we say that the 
control v steers U into U' in t units of time and U' is attainable or reachable from U at time t. 

Definition 1 (Reachable Set): The set of all U' G G attainable from Uq at time t will be denoted 
by R(Uo,t). Also we use the following notation 

K(U ,T) = |J R(U ,t) 

0<t<T 

R(U ) = (J R(U ,t). 

0<t<oo 

We will refer to R(J7q), as the reachable set of Uq. 



Remark 1 From the right invariance of control systems it follows that R(Uq,T) = R(e,T)Uo, 
R(E7 , T ) = R(e,T)U a , and R(U Q ) = R(e)U . Note that R(U ,T) need not be a closed set, we use 
R(Uo,t) to denote its closure. 

Definition 2 (Infamizing Time): Given Up € G, we will define 

t*(U F ) = inf {t > 0| Uf G R(e,t)} 
t*(KU F ) = inf {t > 0| W F G R(e,t), k G K} 

and t* (U) is called the infimizing time. 

From a mathematical point of view, we may identify two goals in this paper: (1) to characterize 
R(e,t) and hence compute t*(U F ), the infimizing time for Uf G G, and (2) to characterize the 
infimizing control sequence v n in (|2|), which in the limit n — ► oo, achieves the transfer time t*(U F ) of 
steering the system (||) from identity e to Uf- From the physics point of view, these results will help 
to establish the minimum time required and the optimal controls (the rf pulse sequence in NMR 
experiments) to achieve desired transfers in a spectroscopy experiment. 
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Figure 1: The panel shows the time optimal path between elements U and V belonging to G. The 
dashed line depicts the fast portion of the path corresponding to movement within the coset KU 
and, in traditional NMR language, corresponds to the pulse and the solid line corresponds to the 
slow portion of the curve connecting different coscts and corresponds to evolution of the couplings. 

3 Time Optimal Control 

The key observation is the following. In the control system (Q), if Up € K then t*{Up) = 0. To 
see this, note that by letting v in (^) be large, we can move on the subgroup K as fast as we wish. In 
the limit as v approaches infinity, we can come arbitrarily close to any point in K in arbitrarily small 
time with almost no effect from the term Xd- By same reasoning for any U € G, t*(U) — t*(kU) for 
k € K. Thus, finding t*(Up) reduces to finding the minimum time to steer the system (0) between 
the cosets Ke and KUf. 

This is illustrated in the Figure 0, where the cosets KU and KV are depicted and the infimizing 
time path between elements U and V belonging to G is shown. The dashed part of the curve 
illustrates the fast motion within the coset. The solid part of the curve corresponds to the drift 
part of the flow ( also known as the evolution of couplings in NMR literature). The minimum 
time problem then corresponds to finding shortest path between these cosets or, in other words, the 
shortest path in the space G/K. 

With this intuitive picture in mind, we now state some lemmas. 

Lemma 1 Let U G G and X : M — > g be a locally bounded measurable function of time. If X n (t) 
converges to X(t) in the sense that 

lim / \\X(t)-X n (t)\\dt = 0, 

then the solution of the differential equation U = X n (t)U at time T converges to the solution of 
U = X(t)U at time T. 

The proof of the above result is a direct consequence of the uniform convergence of the Peano-Baker 
series. We use this to show 
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Lemma 2 For the control system in equation (||), t*(Up) = t*(KUp). 

Proof: We first show that if A; € if, then t*(k) = 0. Because exp{X 1; . . . ,X m }LA = K, given any 
T > there exists v £ C(T), such that the solution of 

m 

D"=E«i(*)-Xi]f. 1/(0) = e 

i=l 

takes on the value k at time T. Now consider the family of control systems 

m 

U=[X d + aJ2v i {at)X i ]U, U(0) = e. 

i=l 

Rescaling time as r = at, we obtain 

rlJJ Y _ m _ 
or a *■ — ' 

i— 1 

Observe that, by Lemma El as a — > oo, C/(r)| r= T = fc or hin^^oo [/(t)| t _ t = fe. Therefore k € 
R(e,T), for all T > 0, implying t*(fc) = 0. 

We now prove the general assertion. Let f*(f/p) = T, we show that if U% = kl/F for k £ K, then 
= T. Because t*(U F ) = T, for any Ti > T,U F ^ R{e,T\), therefore there exists a family 
of control laws v r [0, T{\ such that the corresponding solutions U r (t) to the equation (g) satisfy 
U r (Ti) — > f/p. From the first part of the proof, for any T2 > Ti there exists a control sequence 
v p [Ti, T2] such that the solutions t/ p (£) to the family of control systems 

m 

U = [X d + v?Xi]U, U{Tx) = U F 

i=l 

satisfies U p (t p ) — > C/i, for t p < T2 and i p — > Ti. Using the continuity of the solution of the differential 
equation to its initial condition and Lemma [j], we conclude that there exists a family of control laws 
u™[0, T2] such that the corresponding solutions U n (t) to the family of control systems 

m 

U=[X d + J2 v i X i\ U > U{0) = e 

i=l 

satisfy U n (t n ) -> t/i, for t n < T 2 . Therefore, U x S R(e,T 2 ). Since T 2 > Ti is arbitrary t*(C7 x ) < T x , 
Because Ti > T is also arbitrary, we infer that t*(Ui) < T. This shows that t*(U\) < t*(Up). Now 
reverse the roles of Up and U\ to get the opposite inequality. This proves the claim. Q.E.D. 



Remark 2 The above observation will help us make a bridge between the problem of computing 
t*(Up) and the problem of computing minimum length paths for a related problem which we now 
explain. 



Definition 3 (Adjoint Control System): Let P 6 G. Associated with the control system (0) is 
the right invariant control system 

P = XP, (3) 

where now the control X no longer belongs to the vector space but is restricted to an adjoint orbit 
i.e., X £ Adx(Xd) = {kr l X d k\k G K}. We call such a control system an adjoint control system. 
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For the control system (||), we say that KU f E B(Uo, t') if there exists a control X[0, t'] which steers 
P(0) = Uo to P(t') E KUf in t' units of time. We use the notation 

B(U ,T)= |J B(U ,t). 



0<t<T 



From Lemma we see that B([/o, T) is closed. 
We use 

L*(KU F ) = inf {t > 0| KU F € B(e,t)} 

to denote the minimum time required to steer the system (^) from identity e to the coset KUf- We 
call it the minimum coset time. 



Theorem 2 (Equivalence theorem): The infimizing time t*(Up) for steering the system 

m 

U= [X d + J2 v i X iP 
z=l 

from U(0) — e to [/f is the same as the minimum coset time L*(KUf), for steering the adjoint 
system 

p = xp, x e Ad K {x d ) 

from P(0) = e to i^Z7 F . 

Proof: Let Q E K satisfy the differential equation 

m 

Q=[J2viXi]Q, Q(0) = e. (4) 

Let P E G evolve according to the equation 

P = (Q-^dQ) P, P(0) = e. (5) 

Then observe that 

^p- = [X d + f^v i X i }(QP), Q(0)P(0)=e, 

z— 1 

which is the same evolution equation as that of U, and since U(0) = Q(0)P(0) = e, by the unique- 
ness theorem for the differential equations, U(t) = Q(t)P{t). Therefore, given a solution U{i) of 
equation ( 
equations ( 



) with the initial condition J7(0), there exist unique curves P(t) and Q(t), defined through 
and (|), s atisfying U(t) = Q(t)P(t). Observe that if U{T) = U F then it follows that 
P(T) E KUf- If Uf E R(e,T), then there exists a sequence of control laws v r [0, T] such that the 
corresponding solutions U r (t) of (Q) satisfy U r (T) —* Uf- Therefore, the solutions P r (t) of the 
associated control system (||) satisfy linv^oo P r (T) E KUf- Because B(e,T) is closed, it follows 
that KU F E B(e,T), which implies that L*(KU F ) < t*(U F ). 

To prove the equality observe that if KUf € B(e,T), then there exists a control X[0, T] such that 
the corresponding solution P(t) to (||) satisfies P{T) E KUf- Because X(t) E AdxiXd), we can 
express X{t) as Q(t) XdQ(t). It is well known that we can find a family v r (t) of control laws 
such that the corresponding solution Q r (t) of 

rn 
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satisfy limbec Jq \\Q(t) - Q r (t)\\dt = 0. Hence, lim,^ ^ \\X{t) - (Q r (t))- 1 X d Q r {t)\\dt = 0. 
0^ 



Using Lemma [Lj we claim that the solutions to family of differential equations 

P r = [{Q r r\t)X d Q r {t)]P\ P r (0) = e 

satisfies lim^oo P r (T) 6 KU F . Therefore, t*(KU F ) < T. Since the choice of T was arbitrary, 
it follows t*{KU F ) < L*{KU F ). Because t*{KU F ) = t*(U F ), it follows that t*{U F ) < L*{KU F ). 
Hence the proof. Q.E.D 



Remark 3 The control system (||) evolves on the group G and induces a control system on the 
coset space G/K through the projection map tt. The adjoint control system ([|) is a representation 
of this induced control system. Observe that since ||X|| = f in (||), we can also define L*(KU F ) as 
the infimizing value of < P, P >2 dt for steering the system 

P = 7 IP, 7 > 

from P(0) = e to P(l) E KU F . 



We will now compute t*(U F ) using the properties of the set AdxiXd). Based on the qualitative 
nature of time optimal trajectories of the system (0), we make the following classification. 

f. Riemannian Symmetric Case In addition to Assumption [l], if we have the restriction 
[m, m] C t, then we are in the Riemannian symmetric case as described in the section ^. We 
can further classify this case based on the rank of the symmetric space G/K. 

• Pulse-drift-pulse sequence(characteristic of single-spin systems) In this case, the rank 
of the symmetric space G/K is one. Roughly speaking the trajectories of the infimiz- 
ing control sequence v r (which in the limit r — > 00, achieves the transfer time t*(U F )) 
converge to an impulse (which resembles an impulse of appropriate shape), followed by 
evolution under drift (for time t*(U F )) and a final impulse. 

• Chained Pulse-drift-pulse sequence( characteristic of two-spin system)In this case, 
the rank of the symmetric space G/K is more than one. The trajectories corresponding 
to an infimizing control sequence v r in (^) converge to a chain of " impulse drift impulse" 
pattern. The infimizing time t*(U F ) is the time spent when the system just evolves under 
drift. 



2. Chatter sequence In this case, G/K is no more a Riemannian symmetric case, i.e. [m, m] <f_ i. 
This is a characteristic of more that two-spin systems. 



In this paper we will confine ourselves to the Riemannian symmetric case. The non-symmetric 
case will be treated in detail in a forthcoming paper. 

Pulse-drift-pulse sequence 

We begin with the first case where the rank of the symmetric space G/K is one. It follows from 
Theorem [j] that m = \J a>Q AdxiaX^). In this case, computing t*(U F ) reduces to finding the 
geodesic distance on the homogeneous space G/K. Given the bi-invariant metric <, >g on G, there 
is a corresponding left invariant metric <,>„, on the homogeneous space G/K arising from the 
restriction of <, >g to m. Let L n (j) represent the length of a curve 7 € G/K under the standard 
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induced metric. There is a one-to-one correspondence between the curves {"f(t) € G/K\j(0) = 
0, 7(1) = tt(U f ), L n (<y[0, 1]) = T} and the trajectories of system (|) satisfying {P(0) = e,P(T) € 
JCJ7i?}. Therefore, L*(KU) is the Riemannian distance between n(e) and 7r(t/) under the standard 
metric <,>„. This is computed in the following theorem, which characterizes geodesies on the 
homogeneous space G/K under the standard metric JO]. 

Theorem 3 Let G be a compact Lie group with a bi- invariant metric <,>, and if be a closed 
subgroup. Let g and t denote their Lie algebras with the direct sum decomposition g = m+6, m = t^. 
Consider the right invariant control system 



U= [X d 



^ViXi]U, U EG, U(Q) = e 



where Vi 6 K, € m, and {X^la = {• Suppose G/K is a Riemannian symmetric space of rank 
one, then t*(Up) is the smallest value of a > such that we can solve Uf — Qi ex-p(aXd)Q2 with 

Q u Q 2 eK. 

Proof: By the equivalence theorem t*{Up) = L*(KUp), where L*(KUp) is the minimum time for 
steering the system 

p = xp, x eAd K {x d ) 

from P(0) = e to KUp. Because G/K is a Riemannian symmetric space of rank one, L*(KUf) 
is the Riemannian distance between o and ir(U) under the standard metric <, > n . From J]J], the 
geodesies under the metric <, > n originating from o take the form 7r(exp(rX)) for X £ m. If 
Uf = Qiexp(tXd)Q2 for Q\, Q 2 € K, then tt{Uf) — 7r(exp(t Q^" 1 X £ ;(52)). To see this note 
that Up = (QiQ2)Q2~ 1 ex P(tX d )Q 2 — (Q1Q2) exp(i Q 2 " 1 X (i (52)- Thus, from the above assertion, 
7r(exp(r Q^ 1 X ( iQ 2 )) is the unique geodesic connecting o to tt{Uf) an d has the length L = t. Hence 
the proof. " Q.E.D 



Remark 4 Roughly speaking, the time optimal trajectory (obtained as a limit of the infimizing 
sequence) for the system (|[), which steers the system form [7(0) = e to Uf — Qi exp(aX < i)Q2, takes 
the form e — > Q2 exp(aXc;)(52 Qi exp(aX£i)(52, where the first and last step of this chain takes 
no time, and the time is required for the drift process(second step). 



We now use illustrate these ideas through some examples. 



1 

1 



represent the Pauli 



Corollary 2 Let U e G = SU(2), and let I x = 
spin matrices. Given the unitary evolution of the single-spin system 

U = + vI x ]U, 

where the control v E R. Let U x — exp(—il x t) represent the one-parameter subgroup generated by 
I x . Given Up € SU(2), there exists a unique (3 S [0, 2tt] such that Up = Uiexp[—i(3I x ]U 2 , where 
Ui,U 2 e U x . The infimizing time t*(U F ) = (3. 



Proof: First note that the Lie algebra g = su(2) has the decomposition m = {il y ,il z }, i = {il x }, 
and Adu T (I z ) = m. Observe from corollary [l] that any Uf € SU(2) has a representation Up = 
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Qi exp[— iaI z ]Q2, where Q±, Qi £ U x . Because exp[— itl z ] is periodic with period 47r, the (3 with the 
smallest absolute value for which Up = U\ exp[— i(3I z ]Ui holds, belongs to the interval [— 2ir, 2n]. 
Because —I z £ Adjj (I z ), we can restrict (3 to the interval [0, 2ir]. The proof then follows directly 
from the Theorem |[ Q.E.D. 



" 


-1 


" 




' 








1 








, n z = 








-1 
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Corollary 3 Let O £ G = 50(3), and let fl x = 1 , fl z = -1 represent 
the generators of rotation around x— and z— axis. Consider the control system 

© = [fl z + vfl x }@, 

where the control Let & x = exp(Q, x t) represent the one-parameter subgroup generated by 

fl x . Given 8_f £ SO(3), there exists a unique (3 £ [0, 7r] such that Qp = Qi exp[(3fl x }Q 2 , where 
Qi,Q2 £ ©it- The infimizing time £*(©f) = P- 

Proof: First note that the Lie algebra g = so(3) has the decomposition m = {fl y ,fl z }, t = {£l x }, 
and Ade^^z) = m. Observe that any 6/ £ SO(3) has a representation 6/ = Qi exp[af2 2 ]Q 2 , 
where Q%, Qi £ Q x - Because exp[^O z ] is periodic with period 27r, the proof is on the same lines as 
Corollary |. Q.E.D. 



We now generalize the example to the case where G = SO(n) , the group of n x n orthogonal matrices. 
The Lie algebra is g — So(n), the set of n x n skew-symmetric matrices. The bi- invariant metric on G 
is < f2, il >= tr(il T fl). Consider the following decomposition of g. Let m consists of skew-symmetric 
matrices which are zero except the first row and column and t consists of skew symmetric matrices 
which are zero in the first row and column. Observe that t generates the subgroup SO(n— 1). Then 
we have 



Corollary 4 Let O £ G = SO(n) and let the control system 

m 

e = [n d + Y / v i n l ]e, e(o) = / 

be given, where fid £ m, Vi £ R, fli £ t and fid £ m, such that expjtfid] has period 2ir. Suppose 
that K — expjflij-Lyi = SO(n — 1). Given Op £ SO(n), there exists a unique (3 £ [0, ir] such that 
Of = ©i exp[(3fl x ]Q 2 , where ©i,©2 £ <d x - The infimizing time t*(Q F ) = (3. 

Proof: Observe that Adxifld) = tn and hence the proof is on the same lines as Corollary |[ Q.E.D. 



Chained Pulse-drift-pulse sequence 

Let us now consider the second case in our classification scheme. We now analyze the case when the 
rank of the Ricmannian symmetric space G/K is greater than one. 

Definition: (Weyl Orbit) Given the decomposition g = m + {, let fj C m represent the maximal 
abelian subalgebra containing Xd- We use the notation Ax d = f) f] AdxiXd) to denote the maximal 
commuting set contained in the adjoint orbit of Xd- The set Ax d is called the Weyl orbit of X d - 
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We use c(X d ) = >0,2A = l,^i6 A x d }, to denote the convex hull of the Weyl 

orbit of Xd, with vertices given by the elements of the Weyl orbit of Xd- 

We compute the infimizing time for the system in the following Theorem (^|), which is a 
generalization of the rank one case. 

Remark 5 Recall from corollary (|l|) that, if A = exp(f)), where f) is the maximal abelian subalgebra 
contained in m, then G = KAK. Therefore given any Up E G, we can express Up = Q\ exp(Z)Q 2 = 
QiQ 2 exp(Ad Q2 (Z)), where Q u Q 2 € K and Z € h. Suppose Z = £" =1 Pi x U Pi > 0, X, E A Xd - 
By choosing X(t) to be AdQ 2 (Xi) for /3, units of time we can steer the adjoint control system 
P = X(t)P from the identity to the coset KUp — K ex.p(Ad,Q 2 (Z)). The claim of the following 
theorem is that indeed the fastest way to get to the coset KUp is to flow on the maximal torus, 
Adg 2 (A), Q 2 E K, containing the coset KUp. 

We now state a convexity theorem due to Kostant Jl9| , which is the main idea behind the following 
theorem. 

Theorem 4 (Kostant's Convexity Theorem): Given the direct sum decomposition g = m + t, 
let t) C Tti represent a maximal abelian subalgebra containing Xd E m. Let r : m — > f), be the 
orthogonal projection of m onto rj. Then T : AdxiXd) = c(Xd), where c(Xd) is the convex hull of 
the Weyl orbit of Xd as defined above. 



Theorem 5 (Time Optimal Tori Theorem:) Let G be a compact matrix Lie group and K 
be a closed subgroup with g and 6 their Lie algebras, respectively such that G/K is a Riemannian 
symmetric space. Let the direct sum decomposition g = m + t, such that m = 6 , be given. Consider 
the right invariant control system 

m 

U=[X d + J2 ViXi]U, U eG, U(0) = e, 

i=l 

where Uj E K, Xd E m, {X^la = 6- Then any Up — Q\ exp(aF)(52, where a > 0, Qi, Q 2 E K, and 
Y E c(Xd), belongs to the closure of the reachable set. The infimizing time t*(Up) is the smallest 
value of a > 0, such that we can solve 

Uf = Qi exp(aY)Q 2 , 
where Qi, Q 2 E K and Y belongs to the convex hull t{Xd) . 

We sketch here the outline of a proof. A rigorous proof from a control theoretic viewpoint will 
be presented elsewhere. 

Once again, we will compute t*{Up) by finding the minimum time it takes to steer the system 
P = XP from X(0) — e to the coset KUp. Following the notation of the corollary [l], let A — exp(h) 
denote the maximal torus contained inside G/K. From corollary [l], we know that G/K — AcLk{A). 
Thus it suffices to prove the theorem for Up E A, therefore we will show that if Up E A then T 
is the smallest value of YllLi a i-> a i — such that Up = exp(^™ 1 OjJQ), where Xi E Aj^, then 
t*(Up) = T. It is immediate that in the adjoint control system P = XP, by letting X to be Xi for 
ai units of time we can reach Up in T units of time. All we need to show is that we can reach the 
coset KU p no sooner. Let P(t) be the shortest (or time optimal) trajectory of the adjoint control 
system P = XP that steers P(0) = e to the coset KUp. Let tta ■ G/K — > A, denote the projection 
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on the maximal subgroup, whereby for A\ € A , tta ■ k A\k — > A\ (note that the projection is 
only unique modulo a Weyl group action, hence an explicit choice is required to make the projection 
unique). Let a(t) £ A be a continuous path obtained from the projection of P(t), onto A. Then 
observe the projection tza induces the map it a, '■ P(t) ~ * a(t). Then the evolution of the curve a(t) 
can be expressed as a(t) = £1 a(t), where ft — T(Adf.(X c i)), for some k € K (recall r : m — > f) is 
the orthogonal projection onto fj). Now using Kostant's convexity theorem, we have € c(Xd). 
Therefore we can write a(t) — (J27Li PiXi) a(t), where X{ € Ax d and J27Li Pi = 1 f° r ft > 0- This 
implies that if T — YnLi a i ^ s ^he smallest value for which Uf € A satisfies Uf — ex p(X)^li 
for ai > 0, then the path a(t) will atleast takes T units of time to reach Uf- 

Remark 6 The theorem characterizes B(e,i), the reachable set for the adjoint system. This is 
given by 

KB(e,t) = Kexp(aY)K, 0<a<t 
where Y belongs to the convex hull c(X<i). 



4 Spin Algebra 



The Lie Group which we will be most interested in is SU(2 n ), the special unitary group describing 
the evolution of n coupled spins \. Its Lie algebra su(2") is a 4™ — 1 dimensional space of traceless 
n x n skew-Hermitian matrices. The orthonormal basis which we will use for this space is expressed 
as tensor products of Pauli spin matrices [[l2| 




The matrices (I x , I y , I z ) are the generators for rotation in the two dimensional Hilbert space and 
basis for the Lie algebra of traceless skew-Hermitian matrices su(2). They obey the well known 
commutation relations 

[I x Iy] = il z ; [I y I z ] = il x ; [I z I x ] = il y . 



Then the basis for su(2") takes the form {iB s } where 

n 

B s = 2«- i n(4a) afc % (6) 
fc=l 

a = x,y,or z and 

I ka = 1 <g> • • • <8> I a <8 1 , (7) 

where I a the Pauli matrix appears in the above expression only at the k th position, and 1 the two 
dimensional identity matrix, appears everywhere except at the k th position. cik s is 1 for q of the 
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indices and for the remaining. Note that q > 1 as q — corresponds to the identity matrix and is 
not a part of the algebra. As an example for n = 2 the basis for su(4) takes the form 

5=1 I\x , Ily ,I\z, l2x , l2y , hz 

q = 2 Zllxhx, 2I\ x l2y , 2Ii x l2z 

^hyhxi 2hyl2y,2Ilyl2z 

Zllzhx, Zllzhy, Zllzhz- 

It is important to note that these operators are only normalized for n — 2 as 

tr(B r B s ) = 5 rs 2 n - 2 . 
To fix ideas, lets compute one of these operators explicitly for n = 2 

1 " 
1 

which takes the form 

10 
10 
0-10 
-1 

We will often refer to the algebra of su(2 n ) as the spin algebra. 



5 Optimal Transfer in Two-Spin Systems 

In this section, we will apply our general results on the time optimal control for the specific case of 
a heteronuclear two-spin system. In particular, we consider the following important heteronuclear 
two-spin system discussed in detail in ||. By going to a rotating frame, the free evolution part of 
the Hamiltonian has been reduced to just a scalar coupling evolution. The system then takes the 
following form. 

Let U € ST/ (4), which evolves as 

4 

U = -i{H d + Y,UiHi )U, (8) 

i=l 

where 



H d 


= 2vr J I Z S Z 


H x 


= 2ttI x 


H 2 


= 2nl y 


H 3 


= 2ttS x 


H 4 


= 2irS y , 



hz = 



1 

-1 
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where I x , I y and I z represent operators for the first spin and have the same meaning as I\ x , I± y and 
Iiz, respectively, as explained in previous section [|. Similarly S x , S y , and S z represent operators 
for the second spin and have the same meaning as I 2x , I 2 y and I 2z . The symbol J represents the 
strength of the scalar coupling between I and S. Observe that the subgroup K generated by {Hi}f—i 
is SU(2) x SU(2). 

We first compute the infimizing time for steering the system (^). 

Theorem 6 For the heteronuclear spin system, described by the equation (0), the infimizing time 
t*(Up) is the smallest value of J2%=i a i> a i > 0' sucn that we can solve 

Uf = Qi exp(-i27tJ(aiI x S x + a 2 I y S y + a 3 I z S z ))Q 2 , 

where Qi and Q 2 belong to K. 

Proof: Consider the direct sum decomposition q = m+t, where m = span{I a Sf3\, t = span{I a , Sp}, 
and (a,(3) £ (x,y,z). Then observe [m, m] G t, [m,t] € m, and [t,t] e t. Furthermore, observe that 
^/ 2 s= — {±IzS zl ±I X S X , ±I y Sy}, therefore for a > 0, we have Adxia t{X d ))) = m. Thus the above 
example satisfies all the conditions of the theorem ^. Hence the proof. Q.E.D 

Now we address the question of maximum possible achievable transfer in some given time T. For 
this purpose we define the transfer efficiency. 

Definition 4 (Transfer Efficiency): Given the evolution of the density matrix p(t) = U(t)p(0)W (t), 
where 

U = -i( H d + utHi )U, U(0) = I, 

i=l 

define the transfer efficiency rj(t) to some given target operator F as 

r ] (t) = \\tr(F^U(t)p(0)U\t))\\. 



Remark 7 In the formula for the transfer efficiency, we always assume that the starting operator 
p(0) and the final operator F are both normalized to have norm one (i.e. tr(F'F) = !)■ 



We will now look at the in-phase and anti-phase transfers in the two-spin system, whose evolution 
is given by equation (||). We give here expressions for maximum transfer efficiencies. We first prove 
some lemmas, which will be required in computing transfer efficiencies. 



Lemma 3 let p — 



1 

—i 




and let S be a real diagonal matrix 



oi 
a 2 
a 3 



If ai > aj > afc > 0, where {i,j, k} S {1,2,3} and let U,V E SO(3), then the maximum value of 
llpT/JEV^pH is dj + aj. 
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Proof: Let 



A = 



01 

o 

Jai 



By definition £ = A^A. Using Cauchy Schwartz inequality HpTf/EV^H < IjAUpH ||AC/p||. Observe, 
the maximum value of ||AVp|| is ^/et; + a,j . Therefore Hp^t/SV^I < a.; + a,j. Clearly for appropriate 
choice of U and V, this upper bound is achieved (For example, in case a\ > a 2 > 03, the bound is 
achieved for U and V identity). Hence the result follows. Q.E.D. 



Lemma 4 Consider the function /(ai, 02,0:3) = sin( J7rai) sin( Jita-i) + sin( J7rai) sin( J7ra3). If 
ai, 02, 03 > and 01+02+0:3 = T, where T < jj, then the maximum value of /(oi, 02, 03) is 
2sin(Jira)sin(Jirb), where a + 2b = t and tan(j7ra) = 2tan(J-7r6). 



Proof: Let 

H(ai, 02, 03, A) = sin(j7roi) sin(j7ro2) + sin(j7roi) sin(j7ro3) + A(oi + 02 + 03 — T). 

The necessary condition for optimality gives J^- = , = 0, = 0, which imply respectively 
that 

7r J(cos( J7roi) sin(j7ro2) + cos( J7rai) sin(j7ra3)) + A = (9) 
7rJ(sin(j7roi) cos(j7ro2)) + A = (10) 
7rJ(sin(j7roi) cos( J-7ro3)) + A = (11) 

From equation ([To]) and (p+|), we obtain that either sin(j7roi) = or cos( Jna2)) — cos(J7ra3)). The 
first condition does not give a maxima as it makes / identically zero. The second condition implies 

J7T02 = 2m7T + J7T0 3 . (12) 



Since 02,03 > and 02 + 03 < T < jj, condition (12) is only satisfied for m = 0. Therefore, 
01 = ol2- Now substituting this in (^) and using the equations (JoJ) and ([To|), we get the desired result 
Q.E.D. 



Theorem 7 (Maximum in-phase transfer) Consider the evolution for the heteronuclear IS spin 
system as defined by Equation (||). Let p(0) = x ^ - and F = y . For t < ^j, the maximum 
achievable transfer 

rj*(t) = sin(jTra)sin(jTrb), 
where a + 2b = t and tan(j7ra) = 2tan(J-7r6). For t > the maximum achievable transfer is one. 



Proof: Let 

A(oi, o 2 , o 3 ) = exp(-i27rJ(oi/ 2; S' 2: + a 2 I y S y + a 3 I z S z )). 

From now on we will simply write A(oi, 02, 03) as A. From Theorem^, any unitary propagator Up 
belonging to the set 

3 

R(e,f) = {Q1AQ2I Qi, Q2 e K on > 0, J2 a *- *}j 

i=l 

can be produced by appropriate pulse sequence in (|^). Therefore we will maximize \\tr(F^U(t)p(0)U^ (t))\ 
for U(t) S R(e, i). Let / = exp{il x , il y , il z } and S — exp{iS x ,iS y ,iS z }. By definition, K — S x I. 
In the expression 

r y (t) = ||tr(Q t 1 FtQ 1 AQ 2 p(0)QjAt)||, 
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p(0) commutes with I and F commutes with S, therefore it suffices to restrict Qi and Q2 to / and 
S, respectively. 



Let s denote the subspace spanned by the orthonormal basis {S x , S y , S z } and i denote the subspace 
spanned by the orthonormal basis {I x , I y , I z }. We represent the starting operator p(0) = -^{S x — 

iSy) as a column vector p = -4= [1 — i 0] T in s. The action p(0) — > Q2P(0)Q\ can then be represented 
as p — > where U is a orthogonal matrix. 

Let Pi denote the projection on the subspace i. A simple computation yields that 

Pi (AS x A*') — sin( Jita-i) sin( Jira^)I x 
Pi(ASyA^) = sm(J-Kai) sm(Jna3)I y 
P/(AS' z A t ) = sin(j7ra2) sin(j7ra 3 )/ z . 



We denote the target operator F = ^(I x — il y ) as a column vector ^75 [1 
action p(0) -> P/(AQ 2 p(0)QjAt) can be written as p — > EVp, where 



% 0] T in i. The 



£ = 



sin( Jira2) sin( J7ra3) 







sin( Jirai) sin( J^c^) 







sin( Jirai) sin(j7ra2) 



Therefore we can rewrite r](t) = \\tr(Q\F^Q 1 AQ 2 p(0)Ql At)|j as rj(t) = Wp^UZVpW, where U 
and V are real orthogonal matrices. Using the result of Lemma (^), we get that for sin( Jnai) > 
sin(j7ra2) > sin(j7ra3) > 0, the maximum value of rj(t) is 

sin(j7rai) sin(j7ra2) + sin(j7rai) sin( J^c^) 
2 ' 

Now we maximize the above expression with respect to a±, a 2 , 03 as worked out in Lemma |]to get 
the above result. 

Now we prove the last part of the theorem. Note for t = jj, the maximum achievable transfer 
is one. Because p(0) and F are normalized, this is the maximum possible transfer between these 
operators. If t > jj, say t = T + , we can always arrange matters so that U(T) — e ( by creating 
a propagator U(T/2) — exp(—i2irJ(-jI z S z )) and then creating its inverse exp(i2irJ(^I z S z )) from 
T/2 to T ). In the remaining -jj units of time, we can produce the optimal propagator. 

Q.E.D 

The optimal transfer curve is plotted in comparison with the transfer achieved using the isotropic 
mixing Hamiltonian in the Figure |^. The unitary propagator U (t) in the isotropic mixing Hamilto- 
nian case takes the form 

U(t) = CX pH^V 2 S, + I X S X + IySy)). 



For small mixing times the transfer amplitude achieved by the optimal experiment is up to 12.5 % 
larger than the transfer achieved by isotropic mixing. This is a previously unknown result that will 
find immediate practical applications in NMR spectroscopy. 



Theorem 8 (Maximum anti-phase transfer) Consider the evolution for the heteronuclear IS 
spin system as defined by equation (||). Let p(0) = y/2I z S~ = y/2I z (S x — iS y ) and F = I~ = x ^ y . 

1G 



Optimal Inphase Transfer 




Time in units of 1/J 

Figure 2: The panel shows the comparison between the best achievable transfer (bold curve) and 
the transfer achieved using the isotropic mixing Hamiltonian for the in-phase transfer in 2 spin case. 
On X axis is plotted time in units of 1/J. 

Then, for t < 1/J, the maximum achievable transfer 7]*(t) is 

\\tr(F*U(t)p(p)Ul(t))\\ = sin(Jirt/2). 
For t > j, the maximum achievable transfer is one. 

Proof: Let 

A = cxp(-i27r,J(aiI x S x + a 2 I y S y + a 3 I z S z )). 

From theorem ^| 

U(t) e {Q1AQ2I Qu Q2 e K a, > 0, < t}. 

t=i 

Let S = exp{iS Xl iS y ,iS z } and I = exp{il x ,il y ,il z }. By definition K = S x I . In the expression 
for 

r) = \\tr(FiQ 1 AQ 2 p(0)QiAQ\)\\, 
let Q2 — Q21 x Q2S) where Q21 £ I and Q2S £ S. Let the optimal Q2* = Q21* X Q2S* be such that 

Q2/*p(0)Q 2 /* t = Q2i*IzS~Q 2 i*^ = a z I z S~ + a y I y S~ + a x I x S~ , 
where a x + a y + = 1. Denote 

ry 2 = ||ir(Ft(g 1 AQ 25 )*(V2/ 2 5-)(g^AtQt ) * )|| 
% = !|ir(Ft(Q 1 AQ 2S )*(V2/,5-)(Qt s AtQt)* )|| 
i h = \\tr(F^(Q 1 AQ 2S )*(V2I x S-)(Ql s A^Qlr )||. 
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Then observe that 



V(t) < azVz + a-yVy + a x r\ x . 



We first compute the maximum of r\ z . Let Pj denote the projection on the subspace generated 
by {I x , I y , I z }, then a simple computation yields 



PKA/^At) = 
P / (A/ z 5 z At) = 0. 



sin( Jirai)Iy 
sm{J-Ka2)I x 



Since {I Z S X , I z S y , I Z S Z } forms an orthogonal pair, we can rewrite 

\\tr{F^Q 1 KQ 2S {^I z S-){Ql s ^Q\) )\\ 

i 1 {t) = yu^v P \\, 



as 



where p = [1 — z 0] 5 



i(j7ra 2 ) 
2 








sin( Jxai ) 
2 





and U and V are real orthogonal matrices. From Lemma [|, it follows that the maximum value of 
r\z is 

sin(j7ra2) sin(j7rai) 
2 + 2 ' 

We can compute the maximum of the above expression under the constraint a±+a 2 =t< 1/{J). The 
maximum value of the above expression is obtained for a\ = a 2 . The maximum value is sin(j7rf/2) 
for t < 1/J. Similarly, the maximum value of rj x and r\ y is as above. Since a| + + = 1, we get 
the desired result. 

The final proposition of the theorem has the same proof as in Theorem | Q.E.D. 

The optimal transfer curve for the anti-phase transfer plotted as a function of time measured in 
units of 1/J is shown in the Figure ^.This clearly shows that the transfer efficiency achieved using 
the known mixing sequence [fLOl is optimal. 



6 Conclusion 



In this paper, we presented a mathematical formulation of the problem of finding shortest pulse 
sequences in coherent spectroscopy. We showed how the problem of computing minimum time to 
produce a unitary propagator can be reduced to finding shortest length paths on certain coset spaces. 
A remarkable feature of time optimal control laws is that they are singular, i.e. the control is zero 
most of the time, with impulses in-between. We explicitly computed the shortest transfer times and 
maximum achievable transfer in a given time for the case of heteronuclear two-spin transfers. In a 
forthcoming paper, we plan to extend these results to higher spin systems. 
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